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Abstract 



The fermionic part of the Schrodinger functional of QCD is formulated in 
the lattice regularization with the staggered fermion. The boundary con- 
dition imposed on the staggered fermion field is examined in terms of the 
four-component Dirac spinor. The boundary terms are different from those of 
the Symanzik's theory in the flavor structure due to the species doubling. It 
is argued that, in the case of the homogeneous Dirichlet boundary condition, 
surface divergence does not occur if the link variables of gauge field are intro- 
duced on the original lattice, not on the blocked one. Its application to the 
numerical calculation of the running coupling constant in QCD is discussed. 
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Introduction 



Numerical calculation of the strong coupling constant a s in QCD has been recently 
attempted by several groups. These works involve various physical and technical ingredients, 
namely, charmonium spectrum M , static quark- ant iquark force 0- |1 and the Schrodinger 
functional of SU(3) Yang-Mills theory in finite volume ||- H. To fix the scale, low energy 
physical quantities were measured in the unit of the lattice spacing or in the unit of the 
linear extent of the finite volume. The scale was related to the running coupling constant 
in the (modified) minimal-subtraction scheme by the improved lattice perturbation theory 
H or through the nonperturbatively defined running coupling constant. The basic Monte 
Carlo simulations were all performed in the quenched approximation (the pure Yang-Mills 
theory). As to the charmonium spectrum, however, the potential model has allowed one to 
estimate systematic errors for the approximation 0. 

The desired next step is surely the calculation of a s based on the simulation with dy- 
namical fermion. In this respect, the sea quark effect on the charmonium system has been 
observed by Onogi et al. [[UJ [[11] and their result is consistent with the procedure adopted 
by Elkhadra et al.. Quite recently, a new determination using the T spectrum and includ- 



ing the effect of dynamical quarks was reported by Davies et al. |12[ . These works actually 
suggest that the precise determination of a s through the quarkonium spectrum is promising. 

On the other hand, the finite size scaling technique adopted in the works of Luscher et al. 
can be also applied to QCD. The Schrodinger functional of QCD in the lattice regularization 



has been recently formulated by Sint ||13|| . He examined the boundary condition of fermion 



field in a finite extent of the time direction for the case of the Wilson fermion [ 14| . He 
derived the boundary terms and argued that the resulted boundary terms are not lattice 
artifacts and meaningful at the continuum limit. He also examined the spectrum of the 
Dirac operator of the free fermion and its squared. He has observed that the boundary 
condition forces a minimal frequency in the system and gives a chance to simulate QCD 
with light quarks. 
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Generally speaking, there exists the arbitrariness of choice for the boundary condition 



of the fermion field. Actually, in his original paper on the Schrodinger functional ||15 
Symanzik adopted a different boundary condition from that is derived from the Wilson 
fermion. Concerning the lattice regularization, another type of lattice fermion, the staggered 



fermion [IB|, is known to work. Then a different boundary condition is expected for it. 

As to the staggered fermion, however, the four-fold degeneracy of flavor may cause us 
both good and harm in calculating a s . It may be useful and economical in the simulation 
for the ideal case of four degenerate massless quarks. However, as a final goal, we need to 
put the realistic mass to each quark: a few MeV for the up quark on one side, but about one 
GeV for the charm quark on the other side. The possibility of giving the non-degenerate 
mass term to the staggered fermion has been discussed by several authors |~7|]~ @- But 



it is yet an open question to give such large mass-difference to the staggered fermion in 
actual numerical computations. More technically, it may need a different treatment for the 
staggered fermion of the boundary effect due to the finite extent of the space-time. 

In this paper, we will formulate the fermionic part of the Schrodinger Functional of 
QCD with the staggered fermion through its transfer matrix, which is formulated by 
Thun et al. [|21|, and work out its boundary condition. Then we will discuss the possi- 



bility of getting along with the staggered fermion in the calculation of a s in QCD. 

Transfer Matrix for the Staggered Fermion 

The Schrodinger Functional in lattice QCD can be naturally formulated by the transfer 
matrix as shown by Sint for the case of the Wilson fermion ||13||. The transfer matrix for 



the staggered fermion has been given by H. S. Sharatchandra, H. J. Thun and P. Weisz |2l 
Firstly we will review their results and fix the notation that we use. The lattice spacing a is 
chosen for unity, whereas the lattice sites are labeled by four integers = (m, ri2, n.3, n^) = 
(n, rii). fi stands for a unit vector of the /x-th direction. Since the mass terms are found not 
to be important when one constructs the transfer matrix of the staggered fermion, they are 
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omitted here. Then the action of free massless staggered fermion is given by 

Sf = g ^2vMx{n)\x(n + fi) - x(n - ft)] , 



(1) 



where x( n ) an d x( n ) are one- component Grassmann variables |22|] [23| [21]. Tj^n) is defined 
by 



^(n) = (-l) ni+ - +n "- 1 , rn = l 



(2) 



It is related to the Euclidean Dirac gamma matrices satisfying 



(3) 



by the equation 



T(n)7 M T t (n + p) = r^(n) I 



(4) 



where 



T(n) = 7r i 7 2 n2 7 3 n3 74 n4 • 



(5) 



First of all, we rescale and denote the variables x( n ) an d x( n ) sucri as 



XnA n ) = X{n) , xlAn) = o x{n)V4{n) . 



(6) 



Since rj^n) does not depend on n 4 , we can define 

ik( n ) = Vk{n)r] 4 {n) 



(7) 



and it follows that 



r)' k (n + k) = -rj'kin) 



(8) 



Then the action (|I]) is rewritten as 

S F = J2{xl 4 ( n )Xn 4 +l( n ) + Xn 4 ( n )xl 4 +l(n)} 
n 4 ,n 

+ Y.{xi 4 (nWk( n )Xn 4 ( n + + xl 4 (n + k)r]' k (n) Xn4 (n)} 



(9) 



n 4 ,n k 



The next step is the characteristic one for the staggered fermion: we introduce a set of 
auxiliary fields <p and ^ using the 5 function for Grassmann variable. Namely, we useQ 



where 



n K>«J n %n» - - *«») 

U4 n,n4 



%n 4 N ~ XnS 71 )) = (^ 4 ( n ) - XnS n )) 



Then the first two terms in the action (ft) can be written as 



(10) 



(11) 



EEli ( n )x„ 4+ i {n) + ifl 4 (n)cp n4+1 (n) } . 



(12) 



ri 4 n 



Accordingly, we introduce two independent sets of operators, x, x} an d 0, 0^ satisfying 
the anti-commutation relations as 



{x(n),xV)} = 6 n , 



{0(n) , 0\n')} = S n> 



(13) 



and otherwise zero. On the Fock space spanned by these operators, we can consider the 
coherent states given by 



Xn 4 +1> ( Pn 4 +l 



X.TI4 > Vfl4 



expH{x t (w)Xn 4 +i(w) + 1 (n)<p ri4+1 {n)} |0) 

n 

(0\expJ2{xl 4 (n) X (n) + ifl l4 (n)0(n)} . 



(14) 
(15) 



It can be shown that these coherent states satisfy the following completeness relation. 



n 4 +l 

exp|~Z![Xn 4 (w)x„ 4 +l(w) +^ 4 (n)v? n4+1 (n)]| |Xn 4 +l^n 4 +l) (Xn 4 '^4 



(16) 



We can see that the weight appearing in Eq. (Tl6| ) can be identified with the kinetic term 
in the ^-direction in the form of Eq. flT2]). The remaining terms in the exponential of the 



1 We have adopted the convention for the measure of the Grassmann integral as fdip^dip tpcp^ = 1 



and the abbreviation, 



d<pLd<p n . = U n {dpL (n)d(p n . (n)}. 



action can be arranged into the product of the representatives of the transfer matrices on 
the coherent state basis. The transfer matrix acting on the Fock space is then found to be 

3 

V 

2 

\ it- K=± 

x II : (V H - X(nj) (<f>(n) ~ X f H) : 



f 1 3 

f F = exp -- £ E [x f ( n + kW k (n)^(n) + ^(n^n^n + k) 



xex pf"lEE [f( n + kWMxin) + <£(™K(™)x(™ + *)]}■ (17) 

Tp is hermitian, but not positive definite. The Hamiltonian follows (Tp) 2 and it is positive 
definite. 



Schrodinger Functional 

Hereafter we will consider the staggered fermion in a three-dimensional Euclidean lattice 
with finite volume (2/) 3 and its time development in a finite time interval 2m. The entire 
four dimensional Euclidean lattice can be represented by 

r = { n M = (n, n 4 ) | < n k < 21 (k = 1, 2, 3), 1 < n A < 2m } . (18) 

We assume that the fermionic variables Xi an d ^ obey spatially periodic boundary 
condition. Temporarily, certain configurations in terms of Grassmann variables are assumed 
on initial and final equitime surfaces. Since the Dirac equation is of first order, the boundary 
condition should fix only half of the components of the field variable on each initial or final 
equitime surface. The choice of the half components is intimately related to the way of 
constructing the Schrodinger Functional. 

The wave functional in the Schrodinger representation is naturally introduced as the 
representative of a state vector in the coherent-state basis, 

*[ x t,0t] = / x t^t^\ (19) 



Its dynamics is described by the kernel which acts on an initial wave functional at n± — 1 
and results in a final wave functional at n A = 2m. For the staggered fermion, it is given by 
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^[xLjVLjXIjVi] = (xlm^lm {Tf) 



,t ,„t 



\2m 



Xi>¥>i 



(20) 



It is this kernel which is referred to the Schrodinger functional in the literature. 

The Schrodinger functional can be represented by functional integral as follows. Inserting 
the completeness relation fll6|) in-between the transfer matrices, and integrating over the 
auxiliary fields <p ni+1 , V 9 n 4 ( n 4 = 1> ' ' ' ? 2m — 1) and xL X2m to eliminate the 5- functions, we 
obtain 

, 2m— 1 2m 
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n4=l 



. 2m- 1 



ri4=l 



Xri4 5 ^Ti4 



, 2m— 1 -. 
eX P| - E [Xn 4 ( n )Xn 4 +l( n ) + A( n )^n4+l( n )] } 



/ II K 4 ^n 4 ]exp{-^} 



ri4=2 



(21) 



Sp consists of three parts, 



ntot — q(T) I o i c(°) 



(22) 



If we use the original conjugate variables x(n, n 4 ) for 1 < n 4 < 2m — 1 and x(n, n 4 ) for 
2 < rz-4 < 2m, which are rescaled back from x« 4 ( n ) an d Xn 4 ( n ) by Eq. (H), and we also 
redefine 



Xi(n) = 2774(71)^(71), 



X2m(n) = <fl m {n) 



(23) 



then, these three parts read 

2m-l 4 



n 4 =2 n /i=l 



(24) 



5b -EE 9%( n ? 2m) 2m)[x(n + fc, 2m) - x(n - fc, 2m)] 
™ fe=i z 
1 



S ^(n^n, 2m)x(n, 2m - 1) , 



(25) 



Sfl = E E 1) + *!, 1) - X(« - fc, 1)] 

n fc=l 2 

+ E^( n M n ' 1 W n ' 2 )- 



(26) 
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Boundary Condition in terms of Four-component Spinor 



Starting from the transfer matrix of the staggered fermion, we have obtained the func- 
tional integral form of the Schrodinger Functional for the free quarks. The first part of the 
action, Sp, given by Eq. (E5), is the usual action for the staggered fermion (in finite volume). 



S$ and are boundary terms involving the boundary values of the fermion field, Xi: Xi 
and X2m, X2m, respectively. In order to examine these boundary terms more closely and to 
communicate with the counterparts in the continuum theory, we will express them in terms 
of the four-flavored four-component Dirac spinor [24J] f25| . We follow here the formulation 
given by H. Kluberg-stern et al. 



The original index of the lattice may be written in the form of 

n p = 2x p + , (27) 

where x p is an integer four-vector and p M is a four-vector whose components are either one 
or zero, x^, labels the site on the sublattice with the spacing twice as large as the original 
one. p p points one of the sixteen variables x{2x + p) near to x p making them associated to 
the site on the sublattice. Accordingly, we adopt the notation such as 

X P (x) = x(2x + p), x p (x) = x(2x + p). (28) 

It follows immediately that 

X(2x + p + p) = J2iXp>( x )$P',P+P + x P '(x + p)V,p-a] > ( 29 ) 
p' 

X(2x + p-p) = Y1Xp'{%)$p',p-[l + X P '(x - p)<V,p+/i] , (3°) 
p' 

Then, noting that r q jJi {2x + p) — r]^(p), we can rewrite the action Sf as 

m— 1 

^ = lEEE {Xp(x)( A p)pp' [x P '(x + p) - x P '(x - p)} (31) 
+ X P ( X ) (Al) P(/ [x P > (x + p) + x P > (x-p)- 2 Xp > (x)} } , (32) 

where 



(A/x) PP ' = 77 M (p)[6 p _ A)/ y + <W,p'] , (33) 



( A l) PP ' = V»{p) [<W,p' ~ <W,p'1 • ( 34 ) 

These A^'s can be expressed by the traces over the products of gamma matrices. First 
we introduce matrices 

R P = (RX = \ + P) = \ 7f7 2 P2 73 P3 7f , (35) 



satisfying unitary conditions 



J pp i 



Y,(Rl)a a (R P >) a a = ^{RlR P >} = 5p 

aa 

Ew:«/=ft' (36) 
p 

Then we can show the following relations. 

(A,) pp , = Tr{Rlj,R p ,} , (37) 
(Aj) pp/ = Tr {^75^7574 , (38) 

DW ( a .)pp'(4)^ = (7,)/ ® V = (r,)A a , (39) 

pp' 

£W« a KM R W = (75)/ ® = at) A a • (40) 

pp' 

With this unitary matrices i? , we can transform x P i x ) m t° the four- flavored Dirac spinor 

as 

^a a ^) = ^Lx P {^Rl)a, (41) 



where Greek and Latin suffices can be taken to denote spinor and flavor indices, respectively. 
In terms of the four-component spinor thus defined, we obtainP] 



2 



The normalization in Eq. (41) is determined so as to scale all dimension- full quantities in the 



four component spinor representation by the new lattice spacing b(= 2a) twice as large as a. 
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m— 1 



£4 = 1 X fj. 



(7m®1)V m + -(75®7;7 5 *)A m 



ip(x) 



(42) 



where 



V^tx) = -[^(x + /2) - ip(x - fj,)] , 
A tM ip(x) = ip(x + jl) + ^(x — P) — 2^(x) 



(43) 
(44) 



Now we consider the surface term, , at x± — m = T/b. It can be rewritten in the 
form of 

S ( P = -J2J2^X( P ,o){x,m) (A 4 - A!)( P) o)(pM) x^ijfom - 1) 
35 pp' 



c\X(p, 



x pp' k=l 



(A 



i)V fc + -(Afe) (Pj o)(pM)A fe x (p0) (a;,m) . (45) 



Here we denoted the four-vector p by a set of a three-vector and a definite forth component 
as (p, 0) or (p, 1). It is useful to introduce projectors which act on the space of the four- 
vectors {p} and project those elements with a definite value for the forth component p 4 . We 
define 



(-Po)pp' = <$ pH 



-4,p'+4 > 



(-Pl)pp' — <5p-4,p'-4 



(46) 



The relations below follow the definitions of Eqs. (0) and ([34] 



XX-AaOpct (A/Joy - <Wap'+A + ^p-Ap'-A - <W (A* = 1> ' * * > 4) 



(47) 



E( A mW A pW = <W+A - <W,p'-A (/i = 1, • • • , 4) . (48) 



By using these relations, the above projectors can be expressed by A 4 and A| as follows. 

1 



(-Po)pp' - 2 XXA4WA4 + AJoy , 



(A)«/ = -£( A 4UA 4 -A 



4/crp' 



(49) 
(50) 



Making use of these projectors, Eq. ( [45]) can be written as 
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,(T) 



x p,p' 

1 3 
+ oEEE^( 



P (A 4 -A^)Pi (aj,m-l) 



P (A fe V fc + -A|A fc )P 



X P '{x,m) 



(51) 



With Eqs. 
spinor. 



a: pp' fc=l 

and we can again rewrite it in terms of the four- component 



S ( P = - y £$(x,m)P r 4 il>(x,m-l) 



+ E E ^ m Wr fc v fc + -r|A fc ) j>{x, m) , 

a; fc=l ^ 



(52) 



where Po and Pi are the projectors acting on the spinor and flavor spaces, which corre- 
spond to Po and Pi, respectively: the unitary transformation by R leads to the following 
expressions. 



p = ^r 4 (r 4 + 1*) = i(i <g> l + 7475 ® 7 4 *7 5 *) , 
Pi = ^r 4 (r 4 - it) = l(i ® l - 7475 ® 74 *7 5 *) . 



(53) 
(54) 



Here we have used the anti-commutation relations of T. 



{r^r,,} 
{rj,i*} 



25^ 1®1, 



-25^ 1 <g> 1 , 



0. 



(55) 
(56) 
(57) 



By a similar consideration, can be expressed in the form of 



^ 0) = E^o)Pir 4 ^(*,i) 



+ E E o)Pi(r fc v fc + ir|A fe )v(aj, o) . 

x k=l A 



(58) 



Now we summarize the total action (22) in terms of the four-flavored four-component spinor: 



Sf 1 = ~ E m)P T i ij(x, m - 1) 
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3 ~ 1 

+ E E <fe m)P (r fc V fc + -T 5 k A k )^(x, m) 

x k=l 1 
m—l 4 i 

x 4 =l x M =l z 



+^^( a!) o)TO( a!) i) 



X 

3 



+ £ E ^ 0)P!(r fc V fc + ^A*)^, 0) . (59) 

x k=l 1 

The field components to be fixed by the boundary conditions, xi-> Xi at x 4 = and X2m, X2m 
at x 4 = T/b, can be expressed by the projectors acting on the four- flavored Dirac spinor as 

P 1 ip(x,0) = p (x) , ip(x,0)P 1 = po(x) , 

P ip(x, m) = p T {x) , 4>(x, m)P = p T (x) , (60) 

where pt{x) and pt{x) (t = 0,T) are certain Grassmann valued functions. 

In the classical continuum limit, the coordinate x M and the field ip(x) are scaled by the 
lattice spacing b as 

Vil = bx^ , = pfiil){x) , v(y) = ^p( x ) • ( 61 ) 

Note also that, in the lattice formulation, it is the component Poip(x) at x 4 = +1 that has 
the coupling to po(x), but the same component Poip(x) at x^ = does not exist. In the 
classical continuum limit 6^0, the former comes close to po{x) to have a local coupling. 
We will express this situation with '+0', namely, Po^/(y,+0) = lim-g^PQ-ip^x, +1). The 
same is true for the component P\^/(y) at x 4 = T. Then the action ( |59"|) of the staggered 
fermion reads in the classical continuum limit 



S? = J q T dy 4 d 3 y*(y) [( 7 „ ® 

- | q L A (y, +o)r 4 P! - %(y)r 4 p +o)} 

- | o L d 3 y{r ?T (y)r 4 P 1 ^( 2/ ,T-0)-^(y,T-0)r 4 P 0% (y)} , (62) 



H.B.C. 



where H.B.C. stands for the homogeneous boundary condition with rj(y) and rj(y) vanishing: 
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Pi*(l/,0) = 0, t(y,0)Pi = 0. 

p *(y,T) = o, tf(y,T)Po = 0. (63) 

It is also clear in the lattice formulation that, in the view point of the path integral, the 
integral variables are ^(y) and ^(y) supplemented with the homogeneous boundary condi- 
tion (0), and rj(y) and fj(y) which are localized on the boundary surfaces can be regarded 
as external sources. 

A few remarks are in order here. It is interesting to note that the boundary condition 
for the staggered fermion (|60| ) and the projectors (Q) and (0), are very similar to those 
Symanzik adopted for fermion field in his paper on the Schrodinger representation of the 
renormalizable field theory [fL5|]. According to Symanzik's prescription, the boundary terms 
of the fermion field in a renormalizable theory can be constructed in the Minkowski space 
as follows. We first find out the operator which generates the time reversal transformation 
for fermion. Such an operator may be chosen as 

B = iJ d 3 y^(y) l0 (-z l0 l^(y) = tjd 3 y {%) 7o P_tf(y) - § (y) 7o P + * (y)} , (64) 

where P± = (1 ± i 7 o 7 5 )/2, and the transformation is given by 

i [B,P^(y)] = +P^(y), i[B,P + y(y)] = -P + *(y). (65) 

The component without a change of the sign is referred to "Dirichlet component" and the 
component with a change of the sign is "Neumann component" . First consider imposing the 
homogeneous Dirichlet boundary condition; 

P_^(y,+0) = 0, ^(y,+0)P_ = (66) 

at i/o = 0. (Here the region yo > is referred to "Dirichlet side". ) It can be achieved by 
adding the following boundary term to the action as an interaction. 

B = % J d 3 y +0) 7o P-^(y, -0) - *(y, -0) 7o P+^(y, +0)} . (67) 

To implement the inhomogeneous Dirichlet boundary condition, 
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P-V(y, +0) = V (y) , +0)P_ = ^(y) , (68) 

the Dirichlet component valued functions, rj(y) and fj(y), are introduced as the sources 
which are localized at the boundary surface and the Neumann components in the Dirichlet 
side are made to couple with the sources. Namely, the following source term is added to the 
action: 

B„ = i J d 3 y [*{y, +0) lo P^ V (y) - fj(yft P+*{y, +0)} . (69) 

We can now clearly see the correspondence between the boundary condition for the 
staggered fermion field and the one that is imposed on the continuum fermion field in 
the Symanzik's theory: the time reversal transformation can be associated to the matrix 
± (7475 <8> 7475), and accordingly, the projectors Po,i = (l±747s®7475)/2 enter the boundary 
terms. Except for the structure in the flavor space due to the species doubling, the boundary 
condition is essentially the same. 



Gauge Interaction 

The gauge field part of the Schrodinger Functional of QCD can be formulated just as that 
of the pure Yang-Mills theory given by Luscher et al. || and that of QCD with the Wilson 
fermion given by Sint |13| |. In the case of the staggered fermion x( n )> one can introduce 



the gauge link variables on the lattice of the spacing a. In this case, there is no essential 
difference in constructing the transfer matrix and the Schrodinger Functional from the case 
of the free staggered fermion [ fZTf . 



It is also possible, following Kluberg-stern et al. [55], to define the four- component Dirac 
spinor in a gauge invariant manner. Although the action in terms of the four- component 
Dirac spinor can be obtained only in the expansion with respect to a, this is enough in order 
to read off the boundary condition for the fermion field in the classical continuum limit. The 
gauge invariant four-spinor can be defined by 
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^ a a {x) = -j= J2(R p ) a a U p (x) Xp (x) , (70) 



where 

U p (x) = U 1 (2x) pl U 2 (2x + Pl i) p2 U 3 (2x + Px l + p 2 2) p3 U 4 (2x + Pl l + p 2 2 + p 3 3) p4 . (71) 
The action reads in the classical continuum limit 



H.B.C. 



Sf = J Q dy A jf d 3 y*(y) [( 7/4 ® 1)Z>„] *(y) 

- rf 3 y {f (y, +0)r 4 P 1 r ?0 (y) - r /o (y)r 4 P ^(?/, +0)} 

- | o L d 3 y{r ?T (y)r 4 P 1 vI/( 2/ ,T-0)-^(y,T-0)r 4 P r /T (y)} , (72) 

where 

D^d, + A,(y). (73) 
The boundary condition is the same as Eqs. (|5^). 

Boundary counterterms 

We will next discuss the possible appearance of divergence at the surfaces by interaction. 
According to Symanzik |l5j , since the effect of the boundary surface can be expressed by the 



local surface interaction as Eq. (|67D , it is possible to apply the ordinary prescription for the 
(perturbative) renormalization with the local counterterms, as far as the theory concerning is 
renormalizable in infinite volume. In QCD, if there would appear divergence at the surfaces, 
they could be subtracted by the boundary counterterms which consist of gauge invariant 
local operators of three or less dimensions. For fermions, such terms are only the bilinear 
operators of dimension three and are given in general as 

AS boundary = ( L d 3 y^(y,T)Z T ^(y,T) + t d 3 y $(y, 0)Z o ^(y, 0) , (74) 
Jo Jo 

where Z t (t=0,T) are some matrices. 

To determine Zq and Zt in the lattice regularization with the staggered fermion, we 

should take into account of the symmetry of the staggered fermion in the finite volume 

space. We should consider the following three symmetry transformations: 
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(i) discrete spatial rotation on (i, j)-plain (i, j=l, 2, 3 ) |L9| 
x — > A~ x x ; 

{A^x)i = Xj, (A~ 1 x)j = —Xi, rest unchange, (75) 
ip(x) -> 5(A) <g> P T (A) ^(x) , 

^(x) ->4>(x)S(A)- 1 ® [R T (A))-\ (76) 



with 



= -^(1 - 7< 7i) , R T (A) = ^=( 1 :- 1 *). (77) 



(ii) (74 <S> 7*)-parity transformation JlS 



i/)(x, -> (74 ® 75) ^(-^j ^4) , 

^(cc,x 4 ) -> ^(-aj,x 4 ) (74^)75) • ( 78 ) 



(iii) (75 £g> 7g)-chiral transformation (when massless) 

V>(ar) -> e^ 5 ® 7 ^ ^(x) , 

4>{x) -> ^(z) e^® 7 ^ . (79) 

By these symmetries, the general form of Z is restricted as 

Z = c 4 (7 4 ® 1) + c 5 (7 5 <g> 74 75) + c 4 (74 ® 7*74*) + c 5 (7s ® 7*7s) 

= C (74 ® l)Po + Cl(74 ® l)Pi + C (75 ® 7*75)^0 + Cl( 75 g> 7*7 5 *) Pi . (80) 

where 7 = Sf=i 7i an d c's are arbitrary constants. Then, taking into account of the 
boundary condition (0), the boundary counterterms can be written as 

AS boundary = AS® + AS< T > , (81) 



AS®= / d 3 y{c®r ]o (y)( 1 ^l)P ^(y,+0) + cf^(y,+0)( 1 ^l)P 1 r ]o (y)} 
Jo 

+ [ L d 3 y { <*H(y, +0)( 75 ® 7*75)^(3/, +0) + cS 0) r ?0 (y)( 75 ® 7*7 5 *)A%(2/) } , (82) 
Jo 
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AS^= f d 3 y{c^ ) y(y,T-0)( l4 ®l)P oVT (y) + cPr ]T (y)( l4 ®l)P 1 y(y J T-0)} 
Jo 

+ I d 3 y{^ T) fj T {y)(j 5 ^j*j*)P r] T {y) 
Jo 

+c [ P^(y, T - 0)( 75 ® i*it)P^{y, T - 0) } . (83) 

These expressions give the general form of the boundary counterterms which may appear in 
the lattice regularization with the staggered fermion. 

According to the above expressions, we note that even in the case of the homogeneous 
Dirichlet boundary condition there remain the surface counterterms, which give the coupling 
between the Neumann components. 

AS surface = 4 0) [ L d 3 y *(!/, +0)( 75 ® 7*7 5 *) W?/, +0) 

J 

+ c { P ( L d 3 y^(y,T -0)( 75 <8) 7*7 5 *) A^(y,T - 0) . (84) 
Jo 

The symmetry of the staggered fermion (in terms of the four-component spinor) does not 
seem to be able to exclude these counterterms. This mean that the boundary condition fl63|) 



could not be upheld |L5[ and the Schrodinger Functional of QCD in the lattice regularization 
with staggered fermion would not be well-defined. 

This is not the case, however. To examine this problem, we consider the boundary 
counterterms at finite lattice spacing. Noting P ip(y, +0) = \im-j^Potp(x, +1) and the similar 
equation for Pi^(y,T — 0), we obtain the lattice counterparts of AS SU rface as 

Surface = E +1)(7B ® 7*7 5 * +1) 

x 

+ $(x,m- 1)(7 5 ®7*7;)Pi^(a;,m - 1) 

X 

c (0) 3 

= 4" E E %K 2) X(n, 2) [x(n - fc, 2) - X (n + fc, 2)] 

Z n k=l 

c (T) 3 

H — - — ^2 E ? 7fc( n ' 2 m — 1) x( n ? 2m — l)[x(n — fc, 2m — 1) — x( n + fc? 2m — 1)] . 

^ n k=l 

(85) 

From this expression, we can see the following nature of the surface counterterms. They are 
local and 0(1) in terms of the variables i/j(x) and ip(x) on the lattice with the spacing b. 
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On the other hand, they are derivative terms and 0(a) in terms of the variables x{ n ) an d 
x(n) on the lattice with the spacing a. As a result, if the link variables of gauge field were 
introduced on the lattice with the spacing b, the surface divergence could occur and we 
would need the above counterterms. On the other hand, as far as they are introduced on 
the lattice with the spacing a, these counterterms may not be necessary. 

It is worth noting that the similar thing happens to the self-energy of the staggered 
fermion. As shown in [^l|] |L9| ||20|| , if the link variables are defined on the finer lattice, the 



self-energy correction induced by the gauge interaction is proportional to the bare mass. 
Then the chiral limit is achieved when the bare mass vanishes. However, if they are defined 
on the blocked lattice, a linear-divergent mass term is actually induced and we need a fine 
tuning to obtain massless fermion. 

Note also that the above result is consistent with the case in which the dimensional 
regularization is adopted for the continuum theory with the boundary condition (|60|). In this 
case, a kind of parity becomes symmetry of the theory (including the boundary condition). 
It is given by the following transformation, 

ip{x, x 4 ) -> (74 ® jDipi-x, x 4 ), 

4>(x,x 4 ) ->^(-a;,X4)(74<8>74). (86) 

By virtue of this symmetry, the surface counterterms are eliminated. 

From these considerations, we can expect that as far as the gauge field is introduced on 
the lattice with the spacing a, the boundary counterterms are given in the following form 
in the continuum limit. 

AS bmndary = [ L d 3 y {cj 0) 77 (y)( 7 4 ® l)P *{v, +0) + cf^{y, +0)( 74 ® l)PiVo{v)} 

•> 

+ f L d 3 y {($P*(y, T - 0)( 74 <g> l)P %(l/) + cS T) %(?/)(74 ® l)PMy, T - 0)} . 
Jo 

(87) 
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Quark Masses 



For the staggered fermion, which is equivalent to four-flavored quarks, the degenerate 
mass M can be naturally introduced as 

m i 

S u = E E Mi> a a {x)^{x) = E • (88) 

x 4 =o x „ e r 

Such degenerate mass term for the staggered fermion does not cause any essential mod- 
ification of the above formulation of the Schrodinger Functional. Sm is simply added to 

Qtot 
J F ■ 

Towards the more realistic QCD, we should lift the mass degeneracy of four flavor quarks. 
This possibility for the staggered fermion has been discussed by several authors jL7j]- ||20|| . 
If we adopt the representation for the gamma matrices such as 74 = 03 <8> 1 and 7« = 0"! ® Oi 
for i = 1,2,3, the non-degenerate diagonal mass terms can be expressed as 

M a S ab = ml ab + m A {^l)ab + imn {llll)ab + imvu (74727i)ab , (89) 
In terms of the one-component field x{ x )i h reads 

4 = EEWWCW (90) 

X 

= E {lmx(n) X (n) + ^m 4 (-)" 1+ ^ + "^ 4 (n)x(n) X (n + (-)"4) 
+i\m l2 {-) n ^r ll {n)r l2 {n + (-)" 1 i)x(n) X (n + (-) ni i + R n2 2) 
+i l -m l2i {-) n ^{n)r ]l {n + (-)" 1 i)r ?2 (n + (-) n H + (-)" 2 2) 

x X (n)x(n + (-ri + (-)™ 2 2 + (-)«*)} . (91) 

We can see that the general non-degenerate mass terms lead to the non-local terms extending 
to the next-to-nearest neighbors. They also break the cubic rotational symmetry in the four 
dimensional Euclidean lattice. If we require the spatial cubic symmetry and the locality 
in time (n 4 ) direction, the non- degeneracy is necessarily reduced to two Such mass is 
given by 
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M ab = ml ab + -j=m' (f% 6 , (92) 

where 7 = Sf=i !%■ The two eigenvalues of M ab are given by M = m + m' and m = m — m'. 
In terms of the one- component field x(x), it looks like 

m 

S'm= E E M ^a>)^(*) (93) 

x 4 =0 f 



E (^X(n)x(n) + _L m '(_)n 1+ n 2 +n 3+ n 4 ^-^(^(n^n + (_)«, 



<) ■ (94) 



Note that this mass term, like the degenerate mass (0), does not cause any essential mod- 
ification of the above formulation of the Schrodinger Functional. S'^ is simply added to 

Qtot 

o F ■ 

The degenerate mass term is local and gauge invariant. For the non-degenerate case 
(|94|), the gauge invariance of the mass term can be assured by the suitable insertion of 
the link variables. As shown in |]T9| [p0|1 , if the link variables are defined on the finer 



lattice with the lattice spacing a, the self-energy correction induced by the gauge interaction 
is proportional to the bare mass. Then the chiral limit is achieved when the bare mass 
vanishes. 

Discussion 

In the classical continuum limit of the lattice QCD with the staggered fermion in finite 
volume, we obtained the action (|72|) supplemented by the boundary condition (|3|). We 
also determined the possible form of the boundary counterterms as (|87D in the case that 
the gauge field is introduced on the lattice with the spacing a. This way of introducing the 
gauge field also keeps the chiral property of the staggered fermion |^T[ ||19[ [f20 . 



In the application to the determination of a s in QCD by the finite size scaling technique, 
it seems a practical choice to adopt the homogeneous boundary condition in the numerical 
calculation of a s . For this boundary condition, the boundary counterterms ( j37|) vanish. 



This means that the Schrodinger Functional of QCD with the four-flavored quarks which 
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is supplemented by the homogeneous boundary condition (|63| ) can become finite by the 
renormalization of the gauge coupling constant and the quark mass. As to the quark mass, 
the degenerate mass M will be the first practical choice. The extrapolation to the massless 
limit M —>■ can reduce the number of parameters to be renormalized in the continuum 
limit. Towards the more realistic QCD, we should lift the mass degeneracy of the four-flavor 
quarks. It may be possible to lift the degeneracy between the two "light" quarks and the two 
"heavy" quarks by using the non-degenerate mass term fl94[) . Even though, this is rather 
crude description of the real quark masses. This kind of problem of the nonzero quark mass 
is yet an open question in the practical numerical calculation for QCD in finite volume. 

There are many works remaining to be done. Especially, as the next step, one should 
calculate the contribution of the quark loop to the finite part of the Schrodinger Func- 
tional, Fi[B] in ref. ||, to the first nontrivial order in the weak gauge coupling expansion. 
The calculation in two different regularization schemes would verify the universality of the 
Schrodinger Functional of QCD. 
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